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Abstract 



Radiation Reaction (RR) effects in the interaction of an ultra-intense laser pulse with a thin plasma foil are investigated analytically 
and by two-dimensional (2D3P) Particle-In-Cell (PIC) simulations. It is found that the radiation reaction force leads to a significant 
electron cooling and to an increased spatial bunching of both electrons and ions. A fully relativistic kinetic equation including RR 
effects is discussed and it is shown that RR leads to a contraction of the available phase space volume. The results of our PIC 
. simulations are in qualitative agreement with the predictions of the kinetic theory. 
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1. Introduction 



Q_i Current laser systems may deliver intensities up to 
c/3 '10 22 W cm~ 2 01 and intensities up to 10 26 W crrT 2 are expected 
O jat the Extreme Light Infrastructure (ELI). In such ultrahigh- 
intensity regime and for typical laser wavelength A ~ 0.8 fim 
£^the motion of electrons in the laser field is ultra-relativistic and 
Radiation Reaction (RR) effects may become important. The 
RR force describes the back-action of the radiation emitted by 
an accelerated electron on the electron itself and accounts for 
' jthe loss of the electron energy and momentum due to the emis- 
sion of such radiation. Apart from the need of including RR 
effects in the dynamics of laser-plasma interactions in the ultra- 
\Q relativistic regime, the latter also offers for the first time the 

opportunity to detect RR effects experimentally ]2j,|3|]. 
t— ( ■ In this paper we present an approach to a kinetic descrip- 
tion of laser-plasma interactions where RR effects are included 
via the Landau-Lifshitz (LL) force [4]. Some properties of the 
kinetic equation with RR are discussed and in particular it is 
proved that the RR force leads to a contraction of the phase 
space volume. Then, PIC simulations are used to study RR ef- 
fects on the acceleration of a thin plasma foil in the regime of 
C\3 Radiation Pressure dominance [5]. Numerical simulations 
suggested that Radiation Pressure Acceleration (RPA) becomes 
the dominant mechanism of ion acceleration at intensities ex- 
ceeding 10 23 WcrrT 2 . Such RPA regime is attractive because 
of the foreseen high efficiency, the quasi-monoenergetic fea- 
tures expected in the ion energy spectrum and the possibility 
to achieve a potentially "unlimited" acceleration [0]. Previous 
Particle-In-Cell (PIC) simulations [7] showed signatures of RR 
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effects at intensities exceeding 5 x 10 22 WcirT 2 and increasing 
nonlinearly with the laser intensity. More recent simulations 
studies of RPA both for thick targets [|8j |9[] and ultrathin tar- 
gets [10] suggested that the inclusion of the RR force cools the 
electrons and may improve the quality of the ion spectrum. 

Our approach to the inclusion of RR effects in a PIC code 
has been discussed in detail in Ref.[ 1 1] where one-dimensional 
(ID) simulations of RPA have been also reported. In the present 
paper we report both additional ID simulations and first two- 
dimensional (2D) simulations using parameters similar to those 
of Ref J12ll where, in particular, the impact of a Rayleigh- 
Taylor-like instability on a thin foil acceleration was studied. 

In classical electrodynamics, the effect of RR can be included 
by means of the LL force 
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where v is the electron velocity, y is the relativistic factor, 
r e = e 2 /mc 2 =s 2.8 x 10~ 9 yum is the classical electron radius, 
A = 2jtc/<d is the laser wavelength and we use dimensionless 
quantities as in the PIC code: time, space and momentum are 
normalized in units of co , c<±r l and mc, respectively. Conse- 
quently, EM fields are normalized in units of ma>c/\e\ and den- 
sities in units of the critical density n c = raw 2 /Ane 1 . 

The LL approach holds in the classical framework and quan- 
tum effects are neglected. As pointed out in JTT], the first term 
of the LL force Eq.(Q]i i.e. the one containing the derivatives 
of the electric and magnetic fields, should be neglected because 
its effect is smaller than quantum effects such as the spin force. 
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However, in Sec.© we show the effect of each term of the LL 
force Eq.(Q]i on the rate of change of the phase space volume. 

2. The kinetic equation with Radiation Reaction 

In this section, a fully relativistic kinetic equation that in- 
cludes the RR effects is discussed. We show a few basic prop- 
erties of the kinetic equation pointing out the peculiarities of the 
RR force whose main new feature is that it does not conserve 
the phase-space volume. 

Generalized kinetic equations for non-conservative forces 
and in particular for the RR force have been known since late 
sixties for the Lorentz-Abraham-Dirac (LAD) equation fl3i[l4ll 
and late seventies for the LL equation 1 15]. Recently, the gen- 
eralized kinetic equation with the LL force included has been 
used to study the RR effects on thermal electrons in a magnet- 
ically confined plasma I la) and to develop a set of closed fluid 
equations with RR IU7l4l9ll . In this paper, we give the kinetic 
equation in a non-manifestly covariant form, see EG3] for 
the kinetic equation in a manifestly Lorentz-covariant form. 

The relativistic distribution function / = /(q, p, t) evolves 
according to the collisionless transport equation 



df 

at 



+ V q -(/v) + V p -(/F) = 0, 



(2) 



where q are the spatial coordinates, v = p/y is the three- 
dimensional velocity, y = yl + p 2 is the relativistic factor and 
F = Fl + F« is the mean force due to external and collective 
fields (Fi = -(E + v x B) is the Lorentz force and Fr is given 
in Eq.(Q~|i). Physically, Eq.© implies the conservation of the 
number of particles. 

The new key feature compared to the usual Vlasov equa- 
tion is that for the RR force F« we have V p • F« + 0. Using 
Lagrangian coordinates q(f), p(0> Eq.© can be recast in the 
equivalent form 



d\nf 
dt 



-V p -F. 



(3) 



According to Eq.©, V p • F provides the percentage of variation 
of the distribution function / within the characteristic time scale 
ui~ x . Integrating Eq.Q along its characteristics, we find that the 
distribution function / remains positive as required. 

Introducing the entropy density in the phase space s(q, p, f) = 
-/(q, p, t) In /(q, p, f), from Eq.© we get the equation for the 
evolution of the entropy density 



dt 



+ V q -(sv) + V p -(sF) =/V p -F. 



(4) 



Integrating Eq.© in the phase space, we get the rate of varia- 
tion of the total entropy S (t) 



dS(t) 
dt 



J d 3 qd 3 pfV v -F. 



(5) 



The Lorentz force F^, = -(E + v x B) gives V p • = identi- 
cally thus V p ■ F = V p • Fr. Moreover, the distribution function 
/(q, p, t) is always non-negative / > thus the sign of dS /dt is 
given by V p • F« solely. 



From the LL force Eq.([T} we get j2' 

v*--(?i)(K E -"(^ B -?) 

+ 4y x e) + x B j - 2v • (e X fi) 



E 2 + B 2 
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In a plasma, the kinetic equation is coupled with the Maxwell 
equations for the self-consistent fields 



Pc " e j=e,i J 



(7) 
(8) 



where p c = \e\n c , j c = \e\n c c, J d 3 qd 3 p /y(q, p, f) = Nj is the 
total number of particles for each species (j = e electrons, j - i 
ions) and Zj is the charge of the particle species in units of \e\ 
(for electrons Z e = -1). For a plasma, Eq.© can be recast as 
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The terms of Eq.® proportional to the charge density p and to 
the current density j come from the first term of the LL force 
Eq.dTJ i.e. the term containing the derivatives of the fields. In 
general, these terms can give either a positive or negative con- 
tribution to V p • Fr. The second term of Eq.© i.e. the term 
proportional to (E 2 + B 2 ) has always a negative sign, its ef- 
fect decreases with increasing electron energy and it is typically 
negligible. The third term of Eq.® comes from the strongly 
anisotropic "friction" term of the LL force i.e. the term propor- 
tional to y in Eq.tQ) (see JH for a detailed discussion of this 
term) and dominates in the ultra-relativistic limit y » 1 . 

It is possible to prove l20ll the following statement: for any v 
such that |v| € [0, 1[ then 

E 2 + B 2 

[(v x E) 2 + (v x B) 2 - 2v • (E x B)] + 



2y 2 



> 0,(10) 



therefore according with Eqs.((5l |9]), the terms of the LL force 
Eq.dTJ that do not depend on the derivatives of the fields always 
lead to a contraction of the available phase space volume. In a 
few special cases, the effect of the terms of the LL force Eq.([TJ 
that depend on the derivatives of the fields (i.e. the terms pro- 
portional to p and j in Eq.®) might lead to an expansion of the 
phase space volume. Anyway, their effect should be negligible 
compared to quantum effects as discussed in ll 1 1. 

We show explicitly the contraction of the phase space in the 
special case of a small bunch of electrons interacting with a 
plane wave where collective fields are assumed to be negligi- 
ble compared with the plane wave fields. Assuming an initial 
distribution / = g(q) <5 3 (p - po), from Eqs.(j5] © we have 



dS(t) 
dt 



d 3 qg(q){2 



E 2 + B 2 



y(Po) 

[(v x E) 2 + (v x B) 2 - 2v • (E x B)]} , 



4y(p > 



(11) 



where vo = Po/y(po)- If the electron bunch counter-propagates 
with the plane wave ([vo • (E x B)] < 0) or propagates in the 
transverse direction ([vo ■ (E x B)] = 0), from Eq.dTTb it is clear 
that RR leads to a contraction of the phase space. In particular, 
in the case of counter-propagation (using |E| = |B|, EB = 0) we 
have V p -F s = - (4nr e /3A)4E 2 [2 r (p )|v |(l + |v |) + l/y(Po)]. 
On the other hand, if the bunch propagates in the same direc- 
tion of the plane wave (vo parallel to E x B), then the con- 
tribution of the friction term (proportional to y in Eq.©) be- 
comes comparable with the contribution of the second term 
(proportional to (E 2 + B 2 ) in Eq.®) and we have V p ■ F s = 
-(4nr e /3A) [4E 2 /(1 + |v |) 2 r 3 (po)] which still leads to a con- 
traction of the phase space but with a rate y 4 smaller than the 
case of counter-propagation. This reinforces the evidence of the 
strongly anisotropic features of the LL force Eq.fTJJi (see full for 
further details). 

The physical interpretation of the above properties is that the 
RR force acts as a cooling mechanism for the system: part of 
the energy and momentum are radiated away and the spread in 
both momentum and coordinate space may be reduced. This 
general prediction is confirmed by our PIC simulations (see 
Sec|3) where we found that RR effects lead to both an increased 
bunching in space and to a noticeable cooling of hot electrons. 

Finally, it is worthwhile mentioning that Eq.© is more gen- 
eral than the Vlasov equation but the PIC approach is still valid 
i.e. the PIC approach provides a solution for Eq.© and it not 
limited to the Vlasov equation i20ll . 



3. PIC simulations 

Suitable approximations to the LL force and our approach 
to its inclusion in a PIC code are described in Ref.[ll]. The 
numerical approach is based on the widely used Boris particle 
pusher and it can be implemented in codes of any dimensional- 
ity. Inclusion of RR effects via this method in PIC simulations 
leads to only a ~ 10% increase in CPU time, which may be 
essential to perform large-scale simulations with limited com- 
puting power. 

3.1. ID simulations 

We first report one-dimensional (1D3P) PIC simulations with 
laser and plasma parameters similar to Ref.^. Previous ID 
simulations in this regime have been reported in Ref. lfTlll where 
a detailed comparison with other work is also made. In the 
present paper we review the basic observations in the ID case 
and we include results at intensities higher than those investi- 
gated in Ref.lfTlll. 

The target is a plasma foil of protons with uniform initial 
density «o = 100« f and thickness I — IA where A = 0.8 fim is 
the laser wavelength and T = A/c » 2.67 fs is the laser period. 
In these simulations, the laser pulse front reaches the edge of 
the plasma foil at t — 0, the profile of the laser field amplitude 
has a "trapezoidal" shape in time with one-cycle, sin 2 -function 
rise and fall and a five cycles constant plateau. We considered 
three intensities / = 2.33 x 10 23 W cm" 2 , / = 5.5 x 10 23 W cm -2 



and / = 10 24 WcirT 2 for both Circular (CP) and Linear (LP) 
polarization of the laser pulse. 

In the CP case, we found that RR effects on the ion spectrum 
are negligible even at intensities of / = 10 24 Wcm -2 as shown 
in FigJT] For CP, electrons pile up and the numerical density 
grows exceeding thousand of times the critical density n c . The 
laser pulse does not penetrate deeply into the target (i.e. the 
effective skin depth is a very small fraction of the foil thickness) 
and electrons move in a field much weaker than the vacuum 
field. 

In Ref.j5| it was expected that RR effects in the radiation- 
pressure dominated acceleration of the thin foil would have 
been weak because in this regime the whole foil becomes 
quickly relativistic, hence in the foil frame the laser wavelength 
A' increases and the typical strength of the RR parameter ~ r e /A 
[see Eq.(Q}] decreases. The present case of acceleration with CP 
pulses appears to confirm this picture. The weakness of RR ef- 
fects may also be explained on the basis of the LL equation for 
an electron moving into a plane wave 02111 . As electrons move 
in the forward direction coherently with the foil (while rotating 
in the transverse plane in the CP field) and the amplitude of the 
reflected wave is weak when the foil is strongly relativistic, the 
situation is similar to an electron co-propagating with the plane 
wave at a velocity close to c, for which the LL force almost van- 
ishes II ill . The relativistic motion of the foil also prevents the 
onset of Self-Induced Transparency by increasing the optical 
thickness parameter ( = nno{/n c A in the foil frame (see [22] 
and references therein) . For smaller target thickness, break- 
through of the laser pulse occurs and RR effects are greatly en- 
hanced also for CP ifllll . 

It may be worth noticing that, at the highest intensity consid- 
ered / = 10 24 W cm 2 , in principle one would expect the classi- 
cal approach to RR to break down due to the onset of quantum 
electrodynamics (QED) effects, as discussed in Ref. lll ill . How- 
ever, it can be shown by a direct analysis of the simulation data 
that the threshold condition for QED effect is not violated in the 
CP case. 
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Figure 1 : Ion energy spectra at t = 66T with (top) and without (bottom) RR 
for Circular Polarization. The laser intensity / is 2.33 X 10 23 Wcm~ 2 (yellow), 
5.5 X 10 23 Wcm~ 2 (blue), 10 24 WcirT 2 (red) and the target thickness isf = I A. 
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Figure 2: Ion energy spectra at t = 14T with (top) and without (bottom) RR 
for Linear Polarization. The laser intensity / is 2.33 X 10 23 Wcm~ 2 (yellow), 
5.5 X 10 23 Won" 2 (blue), 10 24 WcnT 2 (red) and the target thickness is t = U. 



For linear polarization (LP), differently from the CP case, we 
found that RR effects are important leading to a reduction of 
the maximum achievable ion energy and to some narrowing of 
the width of the ion spectrum as shown in Figf2] This differ- 
ent dynamics for LP is correlated with the strong longitudinal 
oscillatory motion driven by the oscillating component of the 
j x B force which is suppressed in the CP case. This allows a 
deeper penetration of the laser pulse into the foil with a signifi- 
cant fraction of electrons on the front surface moving in a strong 
electromagnetic field of the same order of vacuum fields fiUl . 
The relative reduction in the ion energy when RR is included is 
close to the percentage of the laser pulse energy which is lost 
as high-energy radiation escaping from the plasma. 

The results for LP (Figf2]i are shown for the same intensity 
values of the CP case (FigQ]) for a direct comparison. How- 
ever, at least for the highest intensity case, the LP results must 
be taken with some caution as the condition for the validity of 
a classical approach tends to be significantly violated. In such 
regime, an analysis based on quantum RR effects might be nec- 
essary 1231I24J]. 

3.2. 2D simulations 

We report preliminary two-dimensional (2D3P) PIC simula- 
tions with laser and plasma parameters similar to Ref . II 1 2TI . To 
the best of our knowledge, this is the first paper reporting results 
of two-dimensional PIC simulations with RR effects included. 

The target is a plasma slab of fully ionized deuterium (Z/A = 
1 /2) of width 40/1, density «o = 169« c and thickness I = 0.5 A. 
The size of the computational box is 95 A x 40A with a spatial 
resolution Ax = Ay = /1/80 and 625 quasi-particles per cell 
corresponding to a total of 8 x 10 7 quasi-particles. The laser 
pulse is s-polarized with the electric field along the z-axis. Its 
normalized amplitude is ao = 320 corresponding to an intensity 
/ = 1 .4 x 1 23 W cirT 2 with a wavelength A = 1 .0 fim and period 
T — A/c « 3.3 fs. The pulse has a Gaussian transverse profile of 
width 20/1 FWHM and a sin 2 longitudinal profile of length 40/1 
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Figure 3: Plots of the 2D PIC simulations at t = 70T. The laser pulse is s- 
polarized with an intensity / = 1.4 X 10 23 Wcm~ 2 and the target thickness is 
t = 0.5/1. From top to bottom, ion «, and electron n e density distributions with 
(left column) and without (right column) RR, longitudinal E x (first row) and 
transverse E z (second row) electric field, ion and electron spectrum with (red) 
and without (blue) RR. 



FWHM. In these simulations, the front of the lase pulse reaches 
the foil at t = 0. 

Comparing the results of our simulations with and without 
RR (see Figf3] we report the results at t — 70T) it is apparent 
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that RR leads to both an increased electron and ion bunching 
and to a strong cooling of electrons. These results are qualita- 
tively consistent with our expectations from the kinetic theory 
that we have discussed in Secf2]and in particular with the pre- 
diction of a contraction of the electrons available phase space 
volume. 

A qualitative understanding of these results can be achieved 
recalling that the RR force Eq.dTJ is mainly a strongly 
anisotropic and non-linear friction-like force that reaches its 
maximum for electrons that counter-propagate with the laser 
pulse IT 1 1 f1 . The backward motion of electrons is thus impeded 
by RR, more electrons and consequently ions are pushed for- 
ward leading to an enhanced clumping that improves the effi- 
ciency of the RPA mechanism. In fact, the ion spectrum with 
RR shows a region between about three hundred and six hun- 
dred Me V with a significant increase in the number of ions com- 
pared to the case without RR (FigO. This picture is confirmed 
by both the enhancement of the longitudinal electric field E x 
and by the formation of denser bunches in the ion density com- 
pared to the case without RR (see Figj3j. However, for linear 
polarization, hot electrons are always generated by the oscillat- 
ing component of the j x B force. The generation of hot elec- 
trons provides a competing acceleration mechanism to RPA and 
ultimately leads to the generation of the fraction of ions with the 
highest energy. The noticeable suppression of the j x B heating 
mechanism due to the RR force therefore leads to a lower maxi- 
mum cut-off energy both in the electron and in the ion spectrum 
(see Fig f3]»- 

These preliminary results for two-dimensional simulations 
with RR effects included suggest that, in the LP case, the trends 
found in one-dimensional simulations hold qualitatively even 
for higher dimensions. More detailed studies and quantitative 
comparisons between one-dimensional and two-dimensional 
PIC simulations are left for forthcoming publications. 

4. Conclusions 

We summarize our results as follows. Radiation Reaction 
effects on the electron dynamics in the interaction of an ultra- 
intense laser pulse with a thin plasma foil were studied analyt- 
ically and by one-dimensional and two-dimensional PIC simu- 
lations. The details of the numerical implementation of the RR 
force in our PIC code were described in Ref. lll ill . 

In one-dimensional simulations, we checked RR effects for 
three different intensities: / = 2.33 x 10 23 WcirT 2 , / = 5.5 X 
10 23 WcirT 2 and / = 10 24 WcirT 2 comparing the results for 
Circular and Linear Polarization of the laser pulse. For CP, 
we found that RR effects are not relevant even at intensity 
of / = 10 24 Wcirr 2 whenever the laser pulse does not break 
through the foil. In contrast, for LP we found that RR effects 
are important reducing the ion energy significantly. 

In two-dimensional simulations, we found that RR reduces 
the j x B heating mechanism leading to a lower maximum cut- 
off energy both in the electron and in the ion spectrum. More- 
over, RR increases the spatial bunching of both electrons and 
ions which are collected into denser clumps compared to the 



case without RR. This might lead to a somewhat beneficial ef- 
fect with a longer and more efficient radiation pressure acceler- 
ation phase whose signature would be an ion energy spectrum 
peaking at an intermediate energy. 

A generalized relativistic kinetic equation including RR ef- 
fects has been discussed and we have shown that RR leads to a 
contraction of the available phase space volume. This predic- 
tion is in qualitative agreement with the results of our PIC sim- 
ulations where we observed both an increased spatial bunching 
and a significant electron cooling as discussed above. 
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